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Abstract 

Using indium as a test case, we investigate the accuracy of the electron- 
phonon coupling calculated with state-of-the-art ab initio and many-body 
theory methods. The ab initio calculations — where electrons are treated in 
the local-density approximation, and phonons and the electron-phonon in- 
teraction are treated within linear response — predict an electron-phonon 
spectral function a 2 F(uj) which translates into a relative tunneling conduc- 
tance that agrees with experiment to within one part in 10 3 . The many-body 
theory calculations — where a 2 F(oj) is extracted from tunneling data by 
means of the McMillan-Rowell tunneling inversion method — provide spec- 
tral functions that depend strongly on details of the inversion process. For the 
the most important moment of a 2 F(uj), the mass-renormalization parameter 
A, we report 0.9±0.1, in contrast to the value 0.805 quoted for nearly three 
decades in the literature. The ab initio calculations also provide the trans- 
port electron-phonon spectral function a 2 r F(oj), from which we calculate the 
resistivity as a function of temperature in good agreement with experiment. 
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I. INTRODUCTION 



Materials with phonon-mediated superconductivity were most intensively studied two to 
three decades ago. In the last decade, especially since the discovery of compounds with 
high transition temperatures, the experimental study of such low-temperature materials has 
steadily declined. In contrast, recent years have seen a steady improvement in computational 
and theoretical methods aimed at describing the electron-phonon coupling in the "old" ma- 
terials. State-of-the-art ab initio methods can now be used to study details of the interaction 
between electrons and phonons and to estimate transition temperatures. The accuracy of 
these calculations, in itself worth investigating, also raises the question of how precisely we 
know the parameters derived from experiments. 

The understanding of phonon-mediated superconductivity relies on a detailed description 
of the coupling between phonons and electrons, most explicit in the electron-phonon spectral 



which phonons scatter electrons on the Fermi surface with an imparted energy u. With 
the addition of an effective Coulomb repulsion, i.e., the Morel- Anderson pseudopotential fi*, 
a 2 F(u) determines all the thermodynamic properties of a phonon-mediated superconductor, 
including the transition temperature To, the critical field, and the specific heat jump at Tc J 
Rather than treating all scattering events equally, one can also weight each scattering event 
according to how much the direction of the electronic velocity changes. This weighting 
results in the transport electron-phonon spectral function a 2 r F{u), which determines the 
transport properties in the normal state. 

First-principles density-functional calculations can be used to study the electronic struc- 
ture, vibrational properties, and electron-phonon coupling in real materials. To calculate 
quantities such as a 2 F{uj) and a 2 r F(u), which involve averages over all phonon modes, the 
density-functional linear-response approach is particularly useful.i^i In this approach, the 
electronic response to atomic displacements is determined self-consistently, and phonon wave 
vectors throughout the Brillouin zone are accessible without having to construct the large su- 
percells needed in finite-difference-based frozen-phonon or generalized supercell methods. 
This approach has been successfully used to study the electron-phonon interaction and su- 




spectral function measures the strength with 
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perconductivity in many simple metals that are suitably treated by the approximations 
inherent in the method, i.e., the local density approximation (LDA) for the electrons and 
the harmonic approximation for the phonons.i'@0 

Alternatively, a 2 F(u) and a 2 r F(u) can be extracted from tunneling experiments and 
optical conductivity data, respectively,0'0 using the Migdal-Eliashberg theory of super- 
conductivity. In particular for a 2 F(u), structure in the tunneling conductance measured 
across metal-insulator-superconductor junctions reflects structure in the superconducting 
gap function A(E) resulting from the interaction of the electrons with the phonons.0ll! 
In the McMillan-Rowell tunneling inversion method,@ the Eliashberg equations are solved 
iteratively to find an a 2 F(u) that is consistent with the measured tunneling spectrum. In 
recent years, improved computational strategies for solving the Eliashberg equations have 
been developed, allowing for more accurate extractions of a 2 F[uj) from tunneling data.@ 

In this contribution, we focus on a 2 F(u) determined from ab initio and many-body 
theory methods to investigate the accuracy with which electron-phonon parameters are 
known. The moments of a 2 F(u), e.g., the electron-phonon mass-renormalization parameter 
A, are quoted for many materials with several digits. Our results show that this is misleading, 
since for indium we find that A can only be given to within approximately 10%. 

Indium is an ideal candidate for our discussion because high-quality tunneling data are 
available. It is possible to fabricate clean tunnel junctions and the electron-phonon coupling 
strength appears large enough to yield a good signal-to-noise ratio.0@ From the theoretical 
point of view, indium also serves as a good test case, since it is a relatively simple metal in 
which relativistic effects are small. The core-valence interaction can be accurately treated 
with a pseudopotential, and since there are only s and p valence electrons, the electronic 
wave functions can be expanded efficiently in plane waves. Furthermore, our results indicate 
that the structural and electronic properties of indium are well described by the LDA and 
that anharmonic effects are small at temperatures near or below the Debye temperature. 

Structure of this paper. The first-principles and many-body theory calculations are out- 
lined in Section II. We refer readers to the references for discussions of the many finer, 
technical points not included here. We present and discuss the results of our calculations in 
Section III, and give concluding remarks in Section IV. 
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II. OUTLINE OF THE CALCULATIONS 



The ab initio procedure consists of three parts: the electronic structure, the vibrational 
properties, and the electron-phonon coupling. 

The electronic structure is calculated in the local density approximation (LDA) of density 
functional theory, by solving the Kohn-Sham equations self-consistently using the Perdew- 
Zunger parameterization of the correlation energy.0 Since it is primarily the valence elec- 
trons that determine the structure and hence interact with the phonons, the core electrons 
are eliminated from the calculation by using a pseudopotential, which is generated by the 
improved Troullier and Martins schemeS The nonlinearity of the exchange and correla- 
tion interaction between the core and valence charge densities is handled with partial core 
correct ions. The Kohn-Sham orbitals are expanded in plane waves with a kinetic energy 
cutoff of 20 Ry. 

Integrations over the Brillouin zone are approximated by sums over discrete sets of It- 
points. The Kohn-Sham orbitals are calculated for 1056 k-points in the irreducible Brillouin 
zone (IBZ). These points, generated with the Monkhorst-Pack scheme JhI originate from a 
mesh of 24 3 k-points in the full Brillouin zone. To accelerate convergence for this metal- 
lic system, we use first-order Hermite-Gaussian smearing with a width of 0.04 Ry.il The 
electronic density of states, in particular the density of states at the Fermi level, N(Ep), is 
calculated more accurately using the linear tetrahedron method.il 

The vibrational properties are determined by calculating the self-consistent first-order 
change in the electron density with respect to atomic displacement s.§ For each phonon wave 
vector q, this change is used to calculate the dynamical matrix, which in turn is diagonalized 
to give the phonon eigenvectors e qi/ and frequencies u>q U [y is a branch label). 

Because the linear-response calculation is the most time-consuming step in the ab initio 
procedure, we calculate the phonons for a relatively small set of 59 q-points in the IBZ. The 
dynamical matrices are obtained on a finer mesh of q-points by a Fourier deconvolution, 
where the calculated dynamical matrices are Fourier-transformed to obtain the real-space 
force constants, which can then be used to form the dynamical matrix at arbitrary q-points. 
The phonon dispersion in the vicinity of q=0 is found to be sensitive to the number of 
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atomic shells included in the force constant model. To ensure the accuracy of the long-range 
force constants, we also do the full linear-response calculation for several small q-points not 
in our original mesh of 59 points. 

The final ab initio step is to calculate the coupling of each phonon to the electron states. 
A phonon qz/ will scatter an electron from a state |nk) to a new state |n'k') with a strength 
determined by the resulting first-order change SV^ CF in the self-consistent potential. For 
atoms of mass M, the electron-phonon matrix elements are given by 



g{nk., n'k', uq) 



h (n'k'|e qi ,-5V q SCF |nk), (1) 



\ 2Muv 

with the restriction k' = k + q. Since only electrons near the Fermi surface can scatter 
via phonons, the average coupling of electrons to a phonon qv is expressed in the doubly- 
constrained Fermi surface average, ((|fiv| 2 ))-^ As with the dynamical matrices, the electron- 
phonon matrices are calculated on the coarse mesh of 59 q-points and then interpolated to 
a denser mesh by means of a Fourier deconvolution. The electron-phonon spectral function, 
which involves coupling to all phonons, is given by@ 

a 2 F(u) =N(E F )J2 Sihcu-hcu^idg^l 2 )). (2) 

In our calculations, the 6 function in Eq. (Q) is replaced by a Gaussian of width 0.5 meV. 

With slight modifications, the above formalism can be used to compute the transport 
spectral function a 2 r F(uj) for the phonon-limited electrical resistivity, p(T), where not all 
scattering events are equally important. For example, forward scattering events do not 
change the direction of the electron velocity v nk and do not contribute to the resistivity. To 
this end an efficiency factor^ 

^nk.n'k' = 1 j n? — , l-Jj 

|v nk r 

is used to weight the electron-phonon matrix elements in the calculations of a 2 r F(u). The 
electrical resistivity is then given byElil 

P{T) - <?N(E F )(v*)2k B TJ hUJ smh 2 (hu;/2k B T) dUJ > (4) 
with Q the cell volume and (v 2 ) the Fermi-surface average of the electron velocity. This is 
a variational solution to the semiclassical Boltzmann equation in which the Fermi surface is 



assumed to undergo a uniform rigid shift in an applied electric field. Here current-current 
vertex corrections are included via the efficiency factor, but only to lowest order. As written, 
Eqs. (^|) and are appropriate for isotropic and nearly isotropic materials since a 2 r F(u) 
is averaged over all directions and (v%) is assumed to be equal to (v 2 )/3. 

We compare the ab initio electron-phonon coupling with experiment both by calculating 
the tunneling conductance from the ab initio a 2 F(uj) and by extracting the experimental 
a 2 F(uj) from tunneling data. The traditional procedure, for which the new computational 
strategies have been developed, is to calculate the experimental a 2 F(u>) from tunneling data 
by solving the Eliashberg gap equations. 

The extraction of a 2 F(u) from tunneling data is done with the McMillan and Rowell 
tunneling inversion method.0 We follow their original prescription: (i) We assume an initial 
value for a 2 F(uj) for which (ii) we adjust fi* to reproduce the experimental superconducting 
gap at zero temperature A , which is 0.541 meV for indium; (iii) we compute the functional 
derivative of the change in the tunneling density of states with respect to a change in the 
assumed a 2 F(u); (iv) we determine the required shift in a 2 F(u) to produce the experimental 
tunneling conductance via a singular-value-decomposition, and (v) we determine the new 
a 2 F(uj) by adding a smoothed shift 5a 2 F{uj) to it. The new a 2 F(u) is then used to begin 
again with step (ii), and the entire process is repeated until it converges. 

While the McMillan-Rowell tunneling inversion procedure is well defined, different com- 
putational strategies can be used to solve the Eliashberg equations. We perform the per- 
turbation theory directly on the imaginary-frequency axis with an energy cutoff of six times 
the maximum phonon frequency, 0J max (beyond which a 2 F(u) = 0), and then perform an 
exact analytic continuation to the real axisJlIizI This method treats the Morel-Anderson 
pseudopotential properly because the sharp cutoff on the imaginary-frequency axis trans- 
lates into a smooth cutoff when analytically continued to the real axis.S In addition, the 
perturbation theory is performed relative to the exact result in the normal state. These de- 
tails are necessary to accurately predict a superconducting Tc from the experimental data. 
The inputs are the experimental tunneling conductance and the superconducting gap at 
zero temperature. The outputs are the electron-phonon spectral function a 2 F(uj) and the 
Morel- Anderson pseudopotential /i*. The transition temperature Tc is then calculated with 
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the T-matrix method of Owen and Scalapinopll rather than using an approximate equation 
such as the McMillan formula. 

The accuracy of the tunneling experiments suffers at low energy, where the signal near 
the gap edge shows a large slope, and at high energy, where the detailed structure in A(E) 
is washed out because it enters the measurement in the form E 2 — A 2 (E). The experi- 
mental data also depends critically on the precise value of the superconducting gap at zero 
temperature, A , because this produces the BCS form for the tunneling conductance; the 
strong-coupling corrections, which are employed to extract a 2 F(u), are the deviations from 
the BCS form. Hence an a 2 F(u) extracted from tunneling data is usually constrained by 
assuming that it has a quadratic dependence at low energy and that it vanishes beyond a 
maximal phonon frequency (some researchers include a quadratic dependence at high en- 
ergies too). Unfortunately, neither the upper limit of the low-energy quadratic behavior, 
nor the exact value of uj max or the frequency dependence near the uj max is known. In prin- 
ciple, uj max should be chosen to be equal to the maximum bulk phonon frequency, but it 
frequently is allowed to be somewhat larger to allow for the effects of interface phonons. 
We adjust the region of quadratic behavior and the maximal phonon frequency in different 
fitting procedures, and impose a linear form on a 2 F(uj) to bring it to zero at uj max . 



III. RESULTS AND DISCUSSION 



The ground-state crystal structure of indium is face-centered tetragonal (fct). As is 
typical with LDA calculations, we find that compared to experiment the equilibrium volume 
is approximately 5% too small. The tetragonal lattice parameters are calculated to be 
a=4.51 A and c=4.84 A; the experimental values are a=4.58 A and c=4.94 A.0 The 
calculated and measured c/ a ratios agree to better than 1%. Our linear-response calculations 
are all performed using the theoretical lattice parameters. 

The calculated electronic density of states (DOS) is plotted in Figure 1. The DOS has 
a free-electron-like behavior at low energies, but develops more structure at higher energies 
where bands cross the Bragg planes. The dashed curve shows the DOS for a free electron gas 
with the same average valence-electron density as indium. The two curves differ significantly, 
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indicating that band-structure and correlation effects within LDA strongly renormalize the 
electron mass in indium. The DOS at Ep is reduced by about 26% compared to the free 
electron value. Unfortunately, it is difficult to compare the calculated DOS directly with 
magnetic susceptibility measurements for indium since these experiments find diamagnetic 
rather than Pauli paramagnetic behavior at low temperatures. 

The inset to Figure 1 shows the DOS near Ep, where the DOS varies only by a few percent 
on the scale of phonon energies. This is in particular true over the range Ep ±6u max , where 
we assume a constant electronic DOS for the tunneling inversion. 

Figure 2 shows the excellent agreement between the measured and calculated phonon 
dispersion curves. The experimental data is taken from neutron diffraction, and is reported 
with an llth-neighbor, 19-parameter Born-von Karman fit.0 Along the direction from Z to 
X for which no direct experimental data are available, we find good agreement between the 
fit to the experiment and our calculated phonon dispersion, though the latter shows more 
structure than the fit. This structure, if real, may be more detailed than can be extracted 
from the available experimental data. 

The electrical resistivity, calculated with Eq. (^), is plotted along with experimental data 
from polycrystalline samples^ in Figure 3. Eq. ([|) is expected to be most accurate in the 
temperature range of about ©_d/5 < T < 20^, with the Debye temperature ®d =129 K.H1 
At very low temperatures anisotropy effects become important, while at high temperatures 
anharmonic effects must be included.il Although the crystal structure of indium is tetrago- 
nal, electrical resistivity measurements on single crystals find nearly the same results along 
the a and c directions. This isotropy also appears in our calculations, where (v 2 ) and (v 2 ) 
differ by less than 5%. 

The dashed curve in Figure 3 is calculated using a 2 F{ui), which is often used as an ap- 
proximation to a 2 r F{uo). The rough agreement between the dashed curve and the measured 
resistivity justifies the approximation in cases when a. 2 r F(u) is not known. However, includ- 
ing the correct weighting with the efficiency factor ^ n k,n'k' brings the calculated resistivity 
into much better agreement with experiment for temperatures up to well above the Debye 
temperature. The transport electron-phonon coupling parameter X tr is found to be 0.74, as 
compared to A = 0.88. The effect of the efficiency factor on the spectral function is shown in 
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the inset of Figure 3, where we also display F(u) scaled to emphasize the strikingly similar 
shape of the phonon density of states and the electron-phonon spectral functions. The ratio 
of a 2 F{u) to F{uj) gradually increases with frequency at about one third of the rate seen in 
leadl 

Figure 4 shows four electron-phonon coupling functions a 2 F{uj): The calculated ab initio 
result and three curves extracted from experimental tunneling data. All three extracted 
curves are based on the same tunneling data taken at T = 0.35 K,0 which should be more 
accurate than the data taken at higher temperatures.El The extracted curves differ in the 
constraints imposed on their low- and high-frequency behavior. The unconstrained curve is 
quadratic for u < 0.5 meV and uses u max = 21 meV. While this curve yields the best fit 
to the relative tunneling conductance, it shows what is believed to be unphysical behavior 
at low and high energies. There is no reason to believe that there is a phonon feature at 
1.5 meV as shown in the unconstrained curve. Rather, that shoulder is most likely an artifact 
related to the accuracy of the voltage (and of Ao) for the experimental data collected at 
low energies. The features at high energy may be real, i.e., structure from either vertex 
corrections or from interface phonons, but most likely they arise from forcing an accurate 
fit to the experimental data at approximately 13 meV above the superconducting gap. We 
estimate that vertex corrections lead to a small reduction of Tc of approximately 0.3%, 
based on a simple integral of ot 2 F(uj) (using the Fermi-surface average C = 0.18).@ This 
result is the same size of effect as seen in lead, so vertex corrections can be safely neglected 
for indium. 

The constrained curve is more strongly restricted in its shape at both low- and high- 
energies: a 2 F(uj) is forced to increase quadratically in u for u < 2 meV and decay linearly 
to zero at the maximum bulk phonon frequency of 16 meV. These constraints eliminate 
what appear to be unphysical features in the unconstrained a 2 F(u), while still fitting the 
tunneling data extremely well. The fourth curve shown in Figure 4 is that of Dynesjll on 
which the three-decade-old value for A is based. His calculation differs from ours in that it 
was performed directly on the real axis, which does not handle fx* properlyil and it was not 
performed relative to the normal state, which would enforce the correct energy dependence 
at high energies. 
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For each of the a 2 F(uj) shown in Figure 4, the differences between the measured and 
calculated tunneling conductance are plotted in Figure 5. It is remarkable that an ab initio 
calculation with only one adjustable parameter (the Morel- Anderson pseudopotential //*, 
adjusted to give the superconducting gap at zero temperature) can fit the experimental 
tunneling conductance to better than one part in 10 3 (the tunneling conductance is on the 
order of 1). Furthermore, we see that the low- and high-frequency features unique to the 
unconstrained curve greatly improve the fit only in the low and high energy ranges. Since 
the experimental data is least accurate in these ranges, it makes sense to constrain the fitting 
procedure to suppress the unphysical features that stem from these ranges. Ideally, the best 
way to proceed would be to use experimental data that has error bars reported with it. 
Such data would allow a maximum entropy technique to be employed to produce the best 
fit a 2 F(uo)B 

Table I describes the curves extracted from tunneling data and the ab initio calcula- 
tion. The electron-phonon spectral functions are characterized by several moments:! (i) the 
electron-phonon mass-renormalization parameter A (twice the first inverse moment), (ii) 
the strength A (the area under the curve), and (iii) the characteristic phonon energy ui n (a 
logarithmic moment). The extracted a 2 F(u) and Morel- Anderson pseudopotential \i* (ad- 
justed to reproduce the experimental A ) are employed to calculate the critical temperature 
Tq with no further adjusting of parameters. All results lie within 5% of the experimental 
Tq- The errors in the tunneling conductance are all rather close (with exception of the 
unconstrained curve), whereas the curves differ significantly in their moments. 

One possible explanation for the wide variation in the moments is the smallness of A 
for indium combined with the experimental uncertainty in the superconducting gap at zero 
temperature, Ao- Since fi* is adjusted to give Ao and materials with small A do not display 
strong features in the tunneling DOS, it is difficult to extract a 2 F(u) to high accuracy. The 
value of a 2 F(uj) in the region between and 3 meV has a large effect on the size of the 
extracted A, but this is the region where the experimental data depends most on the precise 
knowledge of A and the experimental voltage. 

Given the range of A from 0.8 to 1.1 for the different a 2 F(u>), it is not surprising that 
/i* also spans a wide range. Conventional wisdom limits fi* to the range of 0.1 to 0.14 for 
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most materials. In fact, this was a criterion used for choosing junctions in the tunneling 
experiments. Ill However, our ab initio a 2 F(uj) as well as our constrained and unconstrained 
a 2 F{uj) extracted from tunneling data all give \x* larger than the conventional values. Re- 
cent first-principles calculations of \x* suggest that /i* < 0.14 is an artificial limit for some 
simple metals.il Within the standard Eliashberg theory, where a constant electronic DOS 
is assumed, \i* = — fi\n(N(Ep) 6u max )}, and the maximum value is found by letting 
\i become infinite. For indium, this gives a maximum /x* of about 0.25. We get the same 
estimate for a maximal /i* by including the energy dependence of the electronic DOS, i.e., 
the N(E) of Fig. 1, irSi 

with a plasma frequency uj p ~ 12 eV. We expect \i < 1 in indium because it is an s-p metal, 
well described by the free-electron model, so expected values of /i* should be less than 0.2. 
All these estimates indicate that the unconstrained curve is unphysical. 

Even if we discount the unconstrained curve, the remaining values for A differ by up to 
20%. Low-temperature specific heat data can be used to provide an additional estimate for 
A. Using the linear coefficient 7 from experiment^ and our calculated electronic DOS at 
the Fermi level, we estimate A = 0.86, which is close to the ab initio value. This estimate is 
uncertain because the experimental 7 itself is known only to a few percent,^ and furthermore 
the estimate relies on a precise knowledge of the electronic DOS at the Fermi level and on 
the assumption that electron-electron effects do not contribute significantly to the mass 
renormalization. Taken together, these results lead us to conclude that for indium A = 
0.9 ±0.1. 



IV. CONCLUSIONS 

State-of-the-art ab initio methods deliver a very accurate description of the electron- 
phonon coupling in indium: The calculated relative tunneling conductance agrees with ex- 
periment to better than one part in 10 3 ; the calculated intrinsic resistivity as a function of 
temperature is also in good agreement with experiment. The achieved accuracy justifies the 
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approximations invoked: The local density approximation used to calculate the electronic 
structure and the harmonic approximation for the phonons. 

Indium is a good choice for the comparison also because of the high quality experi- 
mental tunneling data. Still, we do not know the strength of the electron-phonon mass- 
renormalization parameter A as well as it would seem from the literature. Based on our 
calculations with state-of-the-art many-body theory methods, we estimate that A can only 
be determined to within 10%, because of uncertainties in the data at low and high energies. 
The uncertainties lead to the question of how to best extract the electron-phonon spectral 
function a 2 F(u>) from experimental data: Is it better to fit the data as precisely as possible 
or to allow for experimental errors at low and high energies by constraining the curve to be 
physically reasonable? All the a 2 F{uj) — ab initio, many-body with and without constraints 
- show the same structure with roughly the same magnitude. The tunneling conductance 
obtained from our ab initio a 2 F(u) is as accurate as the tunneling conductance obtained 
from the most likely spectral function extracted from the experimental data. 

From our study we conclude that the accuracy with which the electron-phonon coupling 
strength is extracted from experiment could be improved. In particular, we hope to motivate 
further experimental work that reports error bars for the tunneling conductance and the 
superconducting gap so that a maximum-entropy technique can be employed to determine 
the best fit a 2 F{u). In materials where vertex corrections are more important, the improved 
accuracy of a 2 F(uo) would also allow the effects of vertex corrections to be observed in the 
multiphonon region. 
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FIGURES 

FIG. 1. Electronic density of states from the ab initio calculation; inset: density of states 
near the Fermi level. The ab initio results (solid curve) are compared to the free-electron density 
of states (dashed curve). The inset shows that the calculated density of states changes by only a 
few percent within the range of Ep ± 6ui max , where uj max is the maximum phonon frequency; this 
justifies the assumption of a constant DOS in this region for the many-body-theory calculations. 

FIG. 2. Phonon dispersion along the high-symmetry directions shown in the inset. The 
experimental data and corresponding fit are from Ref. 30. The ab initio results are calculated on 
a uniform grid of 59 q-points in the irreducible wedge of the BZ and then interpolated to a denser 
mesh by means of a Fourier deconvolution. 

FIG. 3. Electrical resistivity as a function of temperature. The dashed and solid curves 
are calculated using the ab initio results for the electron-phonon spectral function a 2 F(uj) and its 
transport analog a% r F(uj), respectively. The latter agrees significantly better with the experimental 
data.0 The inset shows the strikingly similar shape of the phonon density of states F(u), a 2 F(u>), 
and a 2 r F(u), with the units scaled to display F{uj) with a magnitude similar to that of the spectral 
functions. 

FIG. 4. Electron-phonon spectral functions a 2 F(u>) from ab initio calculations and extracted 
from experimental tunneling results. The three extracted curves are all based on the same tunneling 
data, but differ in the constraints used in the tunneling inversion procedure. The unconstrained 
curve shows spurious behavior at low frequencies and above the maximum phonon frequency. 

FIG. 5. Differences between the measured and calculated tunneling conductances based on the 
spectral functions a 2 F(uj) in Fig. 4. The unconstrained a 2 F(u>) yields errors in the conductance 
that are two orders of magnitude smaller than any of the other a 2 F(oj) curves. 
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TABLES 



error in tunneling conductance 





A 


A (meV) 


Loin (meV) 


* 

M 


To (K) 


max (10 3 ) 


rms (10~ 4 ) 


experiment!! 










3.40 






ab initio 


0.882 


3.00 


5.61 


0.161 


3.31 


1.1 


6 


unconstrained 


1.108 


3.66 


5.20 


0.326 


3.24 


0.04 


0.08 


constrained 


0.984 


3.24 


5.51 


0.224 


3.28 


2.3 


5 


Ref. 17 


0.805 


2.74 


5.84 


0.119 


3.32 


0.9 


7 



TABLE I. Parametric description of the a 2 F(uj) functions from the ab initio calculation and 
extracted from tunneling data. The electron-phonon spectral functions are characterized by sev- 
eral moments: (i) the electron-phonon mass-renormalization parameter A (twice the first inverse 
frequency moment), (ii) the strength A (the area under the curve), and (iii) the characteristic 
phonon energy LOi n (a logarithmic moment). The Morel- Anderson pseudopotential ^* is adjusted 
to reproduce the experimental gap at T=0. Results for the critical temperature Tq, calculated 
from a 2 F(co) and all lie within 5% of the experimental Tc- Also listed are the maximum and 
the RMS errors in the tunneling conductance. 
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Figure 2 — Rudin, Bauer, Liu, Freericks 




Figure 3 — Rudin, Bauer, Liu, Freericks 
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Figure 4 — Rudin, Bauer, Liu, Freericks 
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Figure 5 — Rudin, Bauer, Liu, Freericks 



